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Abstract 

-t— > We present an unconventional approach for addressing the old cosmological constant (CC) 

Q_i problem in a class of F(R, Q) models of modified gravity. For a CC of arbitrary size and sign the 

corresponding cosmological evolution follows an expansion history which strikingly resembles 
that of our real universe. The effects of the large CC are relaxed dynamically and there is 
no fine-tuning at any stage. In addition, the relaxation mechanism alleviates the coincidence 
problem. The upshot is that a large cosmological constant and the observed cosmic expansion 

£■ — 

j- history coexist peacefully in the Relaxed Universe. This model universe can be thought of as 

an interesting preliminary solution of the cosmological constant problem, in the sense that it 
provides a successful dynamical mechanism able to completely avoid the fine-tuning problem 
(the toughest aspect of the CC problem). However, since the Relaxed Universe is formulated 
within the context of modified gravity, it may still suffer of some of the known issues associated 
with these theories, and therefore it can be viewed only as a toy-model proposal towards a final 
solution of the CC problem. 
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1 Introduction 



The triumph of modern physics is largely reflected in our ability to understand the properties of 
the universe and its dynamics from the fundamental physical theories. Yet there is a prediction 
of quantum field theory and string theory which is grossly inconsistent with the existence of the 
observed universe. The vacuum energy density calculated in these theories can be so large in 
absolute value (\p\\ — with Mx around the Planck mass) that the universe would either 

recollapse immediately after the Big Bang or the formation of galaxies that we observe would 
be prevented. This huge gap in our understanding of nature, known as the old cosmological 
constant problem, is one of the most important unresolved issues in fundamental physics [H[2]. The 
cosmological constant (CC) receives numerous contributions from zero point energies of quantum 
fields and phase transitions in the early universe, along with an arbitrary additive constant as 
allowed by general covariance. These effects are of very different orders of magnitude, and all 
of them are many orders of magnitude larger than the value consistent with the observations, 
p\ ps (2.3 x 10 _3 eV) 4 . There is the possibility that these disparate contributions might cancel to 
produce the observed value p\ , but this would entail an unacceptably large fine-tuning. Definitely, 
we have to find another, more natural, solution. 

There have been numerous attempts to solve the old CC problem, e.g. a recent one is [3]. Many, 
if not all of them, and especially those based on dynamical scalar fields, still suffer from severe 
fine-tuning troubles in one form or another p]. Rather than replacing the CC by scalar fields, 
one may think of the CC term as an effective quantity p\ e s evolving with the expansion rate of 
the universe, H. The approach to a dynamical CC in quantum field theory has been recently re- 
emphasized in [1], and it has proven fruitful in tackling the cosmic coincidence problem [J5j. Maybe 
it can also be useful to deal with the old CC problem. Thus, instead of looking for the solution 
in the matter sector, the possibility that it could have a gravitational origin should be seriously 
investigated. This might provide a well behaved "effective CC term" p\ c g(H), viz. one which is 
very large in the early universe (H ~ Mx) but very small at present {H ~ Hq). In this letter, 
we present action functionals of modified gravity which incorporate a very large \p\\ (3> p^) and 
still result in a cosmological expansion with a measurable PAes(H) which is (very) close to the 
dynamics of the benchmark ACDM cosmological model. The striking feature of these modified 
gravity models is that they relax a large CC without fine-tuning. The gravitational mechanism 
for equilibrating the effects of a large CC we refer to as the relaxation mechanism whereas the 
cosmological models in which this mechanism is implemented we call the Relaxed Universe. 

2 Dynamical relaxation of the cosmological term 

An essential ingredient of the new universe is the feedback between the dynamics of space-time 
and the functioning of the relaxation mechanism. The very expansion owing to the existence of 
a large CC sets the relaxation mechanism in motion. The effects of modified gravity equilibrate 
(dynamically) the effects of the large CC. As a result, a realistic universe with a standard sequence 
of radiation dominated, matter dominated and accelerating epochs unfolds automatically. The 
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relaxation mechanism based on the component with an inhomogeneous equation of state (EOS) 
was introduced in [6], while in [7] it was shown that the inhomogeneous EOS is an effective 
description of modified gravity effects. The following step was to construct a viable cosmological 
model. 

In Ref. [8] we have proposed a scenario which achieves this goal in a dynamical way, thereby 
avoiding fine-tuning. In a flat FLRW cosmos, we introduced a dynamical vacuum component, 
which counter-acts the effect of a large \p\\. The total vacuum energy density was found to 
be PAeff = Pa + where f3 is a constant parameter, and the function B representing the 

heart of the relaxation mechanism will follow shortly. First, suppose that p\ ~ — M^ w (with 
Mew = 0(100) GeV) as given by the large (and negative) vacuum contribution from the Higgs 
potential in the electro-weak phase transition. Without counter-measures, the negative p\ would 
let the universe quickly collapse into a catastrophic Big Crunch. To prevent this, one could 
artificially add a fine-tuned CC counter-term. Our relaxation mechanism, though, will do the job 
automatically. Let us illustrate it first with a toy model in the radiation era, where the deceleration 
parameter q = —aa/a 2 is close to but smaller than unity, q < 1. In this epoch, consider (3 > 
and B = H 2 (l — q) given in terms of q and the Hubble rate H = a /a. Clearly, the mechanism 
will stabilize the expansion, because the big p\ < wants to rapidly initiate the collapse, thereby 
pushing the deceleration q closer to 1. However, q can never reach 1 because the second term 
/3/[H 2 (l — q)] > in pAeff will become sufficiently large to compensate the negative CC, finally 
resulting in |pAcff| ^ Ipa| — the CC has been relaxed. And vice versa: if (3/B > would dominate 
over p\ , the value of q would move away from q ~ 1 towards smaller values due to the large positive 
vacuum energy. Consequently, (3/[H 2 {l — q)] will decrease and q will reach a stable point where 
PA and (3/B compensate each other. The resulting universe will expand like a radiation cosmos 
since small changes of q around q ~ 1 are sufficient to keep both terms of PAeff i n equilibrium. 

3 A more realistic model 

To apply the CC relaxation mechanism to the whole Big Bang expansion history, consider the 
more sophisticated B-function in pAeff above 

B= 2 -R 2 + \g + {yR)\ (1) 

involving only the Ricci scalar R = &H 2 {1 — q) and the GauB-Bonnet term Q = —24H 4 q, both 
computed in the FLRW metric. Thus, B reads explicitly 

B = 2AH 4 (q - i \{q - 2) + H e [6 y{\ - q)f . (2) 

It is easy to see that the second term, proportional to H G (1 — q) s , will be responsible for the 
CC relaxation in the radiation regime (q — 1), where the Hubble rate was large. Since, however, 
H decreases with time, the first term oc H 4 (q — 1/2), with a lower power of H, will eventually 
dominate and relax the CC in the matter era, where q ~ 1/2. The arguments for the dynamical 
stabilization in pAeff are exactly the same as in the toy model above. Another novelty is the fixed 
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parameter y ~ H e q^ 3 , which is related to the Hubble rate -ff e q ~ 10 5 Hq at the time of matter- 
radiation equality, where Hq ~ 10~ 42 GeV is the current value. In the late accelerating universe, 
q < and the tiny expansion rate H ~ Hq makes B ~ if 4 the relevant term. In summary, the 
smallness of the denominator B in /?Aeff during all cosmological epochs is the clue to the relaxation 
mechanism. An additional physical argument in favor of it is the growth of inhomogeneities, which 
is in accord with observations [9]. 

This is a good place to demonstrate the absence of fine-tuning in our mechanism. Let us 
ignore 0(1) factors in the Friedmann equation p c = Pm + Pk + fi/H A at late times. Since pa and 
the compensating term (3/H are much larger than the matter density p m and the critical energy 
density p c , we find, in very good approximation, H 4 = \(3/pa\ at late times. Obviously, small 
changes in the parameter (3 will induce only small changes in the Hubble rate H. Hence, it suffices 
to fix the order of magnitude of f3 to obtain a small H ~ Hq. This has nothing to do with fine 
tuning. For instance, take a typical GUT value |pa| — M^, with Mx = 10 16 GeV; then (3 = M 8 
with M ~ 10 -4 eV (of order of a light neutrino mass) would do. An essential result thus follows: 
large \p\\ guarantees small asymptotic H ~ Hq for reasonable M. 

4 An action functional formulation of the relaxation mechanism 

Next we implement the CC relaxation mechanism, here for the first time in the modified gravity 
setup and in the metric formalism [10\ 111] . The crucial part of the model is a function F(R,Q) 
of the Ricci scalar R and the GauB-Bonnet invariant Q. The complete action functional of our 
cosmological model is the following: 
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PA- /3F(R,g)+C mat 



(3) 



where £ m at is the matter Lagrangian, and pa an arbitrarily large cosmological constant. Gn is 
the Newton constant and the parameter f3 will also be dimensionful in general. The variational 
principle 5S/5g ab = then leads straightforwardly to the Einstein equations, 



Gab = SuGn [gab PA + ^PEab + T a b] , 



(4) 



involving the Einstein tensor G a b = Rab — \9abR an d the energy-momentum tensor T a b of matter. 
Additionally, there is a new tensor E a b coming solely from the F(R, Q) term in the action ([3]). On 
a spatially flat FLRW background with line element ds 2 = dt 2 — a 2 (t)dx 2 and scale factor a(i), 
the tensor components in Eq. (Uh are given by G° = —3H 2 , G l j = —5 l j(2H + 3H 2 ) and 



F(R, Q) - 3(H + H 2 )F R + 



- 12H 2 (H + H 2 )F G + 12H 3 F g 

-F(R, Q)-{H + 3H 2 )F R + 2HF R + F R 



+ (H + H 2 )(8HF G - \2H 2 F G ) + 4H 2 F g 



(5) 



(6) 
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where F R and are partial derivatives of F with respect to R and G, respectively, and H = 
-H 2 (q+l). 

The matter sector is described by the standard energy-momentum tensor of an ideal fluid, with 
proper energy density p = T° and pressure p = —T\/3 (summed repeated indices). Similarly, the 
F(R, ^-functional induces the effective energy density and pressure 

p F = 2(3E° and p F = -2 (3 E\/3 . (7) 

Since G a b and E a b are covariantly conserved, T a b is conserved, too. Therefore, the Bianchi identity 
on the FLRW background, p + 2>H(pA-p) = 0, is valid for both matter (p m ,p m ) and F-components 
(pf,Pf), and hence also for (pA e s, PAeft), where p\ e g = p\ + pp is the total energy density of the 
effective vacuum sector and pAeS = — Pa +Pf denotes the corresponding pressure. The generalized 
Friedmann-Lemaitre equations read: 

3H 2 = 87TG N (pm + Pr + PAes), (8) 
H 2 (q- 1/2) = 4 7 rG iV (Pr+PAeff). (9) 

They contain the energy densities of pressureless matter p m = p^a -3 and radiation p r = p®a~ A 
with pressure p r = p r /3, and the effective vacuum components. 

Motivated by the toy- model discussion above, the function F in our setup is taken in the 
simple form F(R, Q) = 1/B, with B given in Eq. Q. The key to the relaxation mechanism in this 
functional framework is the observation that every derivative of F in the expressions ^ and ([§ 
yields another factor £? . Thus the induced quantities pf and pf adopt, in general, the structure 

Ni(H,q,q) N 2 (H,q,q,q) 
PF = ^3 , PF = ^ , (10) 

where Ni% are functions not proportional to B. Correspondingly, the effective vacuum energy 
density adopts the form pAeff = PA + Ni /B 3 and, in complete analogy with the toy- model discus- 
sion, it can be small because the large p\ can be dynamically compensated by letting B become 
sufficiently small (but non-zero), i.e. B — > - hereafter referred to as "relaxation condition". 

Note that, at the equilibrium point (e.g. q ~ 1 in the radiation era), both terms in pAeff 
are almost equal to each other, apart from opposite signs. Therefore, pf behaves approximately 
as another (large) cosmological constant. Obviously, PAeff is not constant in general, because 
Pf = PF(H,q,q) is time-dependent. In addition, there is a corresponding compensation of the 
pressure terms pf and —pa in PAcff- Finally, since the functions §5§ and Q differ from each other, 
the effective EOS, to = PAes/pAeS, will be in general a non-trivial function of time or redshift. This 
should have phenomenological implications. 

5 Analytical and numerical analysis of the model 

In the following, we will describe approximately the cosmic evolution by making use of the CC 
relaxation condition, which is equivalent to |pAcfr| *C \pa\- Remember that B never vanishes. Our 
analysis will be supported by exact numerical results (cf. Fig. [TJ obtained by solving directly the 
Friedmann equation ^ with pp given in Eq. ([7]). 
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Let us start in the matter era where q ~ 1/2 and H 2 ~ p m oc a 3 . Applying the relaxation 
condition in Eq. Q leads to (q — 1/2) oc H 2 , and Eq. (KM) renders 



H 2 (q - 1/2) = ^G N { Pr + p Aeff ) oc i? 4 oc a" 6 , (11) 

whereupon p\ e S = ~Pr/3 + c\ a -6 , with ci a constant; and then via the Bianchi identity for 
{pAeS, PAcs) we find pAeff = C2 a~ 3 — p r + c\ a -6 , the effective vacuum energy density, with c 2 
another constant. Thence 

_ PAeff _ -pr/3 + cia^ 

UJ 3^ l^J 

PAeff c 2 a d - p r + a a 
is the corresponding EOS, which interpolates between dust matter (uj — > 0) at late times and 
radiation (a; — > 1/3) in the early matter era. Depending on the integration constant c 2 a pole 
might occur in uj when pAefT changes its sign, see Fig. [T} It is worth noticing that pAeff behaves like 
dark matter in this epoch, and we could speculate on incorporating both dark matter and dark 
energy into PAeff- Apart from that, we remark the approximate tracking relation PAeff oc p m , which 
can be considered a cornerstone for solving the coincidence problem. 

We expect that the relaxation mechanism became active only as of the early radiation era, 
because F = 1/B ~ in Eq. ^ was sub-dominant during the preceding periods of inflation 
and reheating, where q was far from 1. However, deep in the radiation epoch, before the matter 
era, we find q ~ 1 and F becomes very important. Then H 2 ~ p r oc a -4 implies R = &H 2 {\ — q) oc 

4 

Hi oc a 
relation 



H 3 oc a s upon using the relaxation condition in Eq. (|2l). Moreover, Eqs. (JSh and (|9l) lead to the 



3H 2 {q - 1) = AttG n [3p Aeff - (p AcS + p m )] oc a" 8 / 3 , (13) 

and thus to p\eS = (Pm + PAcff)/3 + csa~ 8 ^ 3 , for constant C3. From here the Bianchi identity 
determines p\ e s = C4 a -4 — p m — (9/4) C3 a~ 8 / 3 , with a new a constant C4. Accordingly, the EOS 
interpolates between radiation at early times and dust matter at later times: 

= PAeS = 1 1 + (3c 3 /4c 4 ) a 4 / 3 

PAcff 3 1 - {p° m /c 4 ) a - (9c 3 /4c 4 ) a 4 / 3 ' 

Again a pole in oj (visible in Fig. [TJ can appear depending on C4. As in the case of the matter epoch, 
it entails no physical singularity since all the energy densities remain finite; it merely reflects the 
change of sign of pAeff- Remarkably enough, we encounter the tracking property PAeff oc Pr oc a~ 4 
also in the radiation era. 

At late times, the universe leaves the matter epoch (p m ~ 1/a 3 — > 0) and gradually enters the 
"dark energy (DE) era". As q departs now significantly from 1/2 or 1, the last dynamical resource 
left to the relaxation mechanism for compensating the big p\, is to choose a very low value of H, 
whereby the DE epoch takes over somewhere near our time (cf. Fig. [I]) - "cosmic coincidence" . As 
the lowest power of H becomes dominant in Eq. ([2]), we have B ~ 2AH 4 (q — 1/2) (q — 2) and the 
full expression for pp in this regime can be computed in closed form: 

PF _ 5q 2 -3q-5 q (4q 2 - lQg + 7) 

P 2m{2q 2 - 5q + 2) 2 2H 5 (2q 2 - 5q + 2) 3 ' 1 J 
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2 3 
log 10 (z+l) 



Figure 1: Deceleration parameter q, effective EOS u) = PAes/pAeS, and relative energy densi- 
ties O ra = Pn/pc of the effective vacuum pA c s (orange thick curve), matter p m (black dashed- 
dotted) and radiation p r (red dashed) as functions of the redshift z. The thick orange curve in the 
q plot corresponds to the relaxation model, and the black dashed-dotted curve to ACDM. Model 
F = 1/B, with inputs y = 7 x 1(T 4 F ~ 2/3 , p A = -W 60 p° c , 0^ = 0.27, = 10" 4 , q rj -0.6, 
q = -0.5 H Q , with p° c = 3H$/(8irG N ). 
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This equation holds for the recent, current and asymptotic future times, i.e. whenever the universe 
is DE or matter dominated. With H — > 0, the previous expression would grow without limit, unless 
the numerators eventually tend to zero, too. Thus, the accelerating asymptotic solution ensues for 
q — > and q — > (3 — vT09)/10 —0.74, see Fig. [I] There are actually three solutions for t — > oo: 
a phantom phase, a de Sitter phase with constant H and an asymptotic power-law expansion 
a(t) oc t r with r > 1. The previously found solution is of this sort, with r — > l/(q + 1) ~ 3.91. 

It is not difficult to work out generalized relaxation models, e.g. of the form F = R n / B m 
with m,n > 0. The crucial ingredient is the function B in the denominator of F. Generally, this 
yields induced terms pf,Pf oc B~ s with s > as a result of Eqs. 
introduce more factors of B in the denominators. One can readily show that the CC relaxation and 
tracking properties in the matter and radiation eras follow again from the relaxation condition. 
Furthermore, the parameter f3 is related to a power of a mass scale M as \f3\ = M 4 ~ 2n+4m . The 
size of M can be estimated by applying the approximations H,q ~ Hq, q G [—1,0], together 
with the condition \p\\ ~ \pp\ ~ \P F\, at H ~ Ho. Again we take \p\\ ~ M x for the standard 
GUT energy density (Mx = 10 16 GeV). In the case (n,m) = (0,1), we recover the previous 
result M 8 = |/3| = 0{p A H^), i.e. M ~ 10- 4 eV. Similarly, we obtain M ~ 100 MeV~ A QC d, 
for (n, m) = (3,2); and M ~ Mx, for (n,m) = (2,1). Remarkably, the mass parameter M of 
the relaxation mechanism is not only completely free from fine-tuning problems, it also lies in a 
perfectly reasonable range of particle physics masses, possibly related to neutrinos, QCD or even 
GUT models. 

6 Conclusions 

To summarize, in this letter, we have unveiled a whole class of modified gravity action models 
which, despite holding an arbitrarily large CC at all times since the early epochs, display a very 
small "effective CC" at present, and without ever needing fine tuning. Like many modified gravity 
models our approach is not complete because of the existence of extra degrees of freedom and 
related problems, which require further considerations, see e.g. I12j. Nevertheless, the Relaxed 
Universe scenario is not just another example for inducing late time cosmic acceleration, for the 
gravity modifications needed to dynamically compensate the large CC are crucial during the entire 
cosmic history. The mechanism also predicts remarkable tracking properties, and characteristic 
dynamical features in our recent past, which alleviate the coincidence problem. A more detailed 
account of this framework will be presented elsewhere. 
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